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ABSTRACT 


An analysis is done of the motion of an tmderwater cradle for submerged 
recovery of the DSRV Alvin. 

The cradle is modeled as a three degree of freedom system with four cor- 
ner springs and includes heave pitch and roll motions. The surface ship Lulu 
is taken as the driving force that couples energy into the cradle through the 
springs. Hydrodynamic damping values are found by a 64 section finite ele- 
ment analysis and are compared to corresponding linear approximations. 

Solution of the resulting equations was accomplished using a Runge-Kutta 
numerical integration methcd developed at M.I.T. 

The SOlutiOn Showed that four sort corner springs in the range soi 20-5¢ 
eevee substantially decrease the transmicced ship motion. Additionally, 
the results indicate that the pitch and roll natural frequency is higher than 
the heave natural frequency and consequently is the limiting design considera- 
ron. 

The results of the analysis are discussed and several possible approaches 
are considered for producing a soft spring constant while at the same time 


supporting the high static load of the cradle weight. 
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INTRODUCTION 


The handling of heavy objects through the air-sea interface is not a 
new or unique problem and has been studied many times in the past. The 
deep submergence research vehicle (DSRV) Alvin owned by the U.S. Navy and 
operated by Woods Hole Oceanographic Institute (WHOI) euros t likely the 
DSRV with the world's best record for successful dives completed. The 
submersible is handled by its support ship, a catamaran named Lulu and is 
transported to and from the open ocean dive site between Lulu's twin hulls. 

The launching and retrieving operation is accomplished with a great 
deal of care, Alvin being initially positioned on her 30-ton capacity cradle 
ee aft of amidships. The cradle, (Figure 1), is supported at its four 
corners by chain hoists which are two-blocked such that the cradle 1s at 
deck level during transit. This allows easy maintenance, repairs and man- 
mine of Alvin. | 

During launching, the catamaran lies with its stern to the seas, the 
cradle is then lowered and the floating submersible clears Lulu with the aid 
of Alvin's screw and numerous snubber lines. Recovery is the reverse of the 
operation, Alvin surfacing well away from Lulu, closing Lulu, and finally 
driving between Lulu's hulls and being manually positioned over the sub- 
merged cradle. Last, when the wave timing is correct, the cradle- 
Submersible combination is raised clear of the water. 

This method of launching and retrieval has allowed reasonably routine 
operations to occur in sea states up to three, with occasional sea state 


four retrieval due to emergency conditions. 
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Figure l. 
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Perot e-cin, 1ntc a discussien of Sea state probabilities, it is well . 
frown that along the north Atlantic off Cape Cod the sea state will frequent- 
ly go from state three to state four due to a prevailing southwest wind which 
increases from ten knots in the morning to about 20 in the afternoon. Sea 
state four makes for a very difficult recovery should it occur during the 
time of a dive. The problem here really is that once the Alvin is in posi- 
tion to be hoisted SSeeee the relative motion between Alvin, the cradle, and 
Lulu could cause very high and damaging shock loading on Alvin, the hoist, 
or the cradle. This loading could possibly result in severe damage to Alvin, 
or even another loss. 

Clearly then, a smooth docking to a stationary cradle followed by a 
rapid hoisting through the air-sea interface would be highly desirable and 
would virtuallv eliminate shock loading and consequent ddifiage. 

Fortunately, the Alvin's cradle consists of four 120' lengths of chain 
and four synchronized wildeats. If the cradle were lowered to approximately 
100' below the surface, the surface wave effect would be virtually inefec- 
tive on Alvin and it would be possible to dock or depart from the cradle 
without the fear of meen cur face concitions. The problem which appears to 
be the limiting factor is the motion of the cradle itself. This motion is 
nearly the same as that of Lulu due as the inelastic chains which connect 
the two. Some work has been done to show the transmitted surface motion 
that the cradle would feel. It has been generally concluded that if members 
more elastic than chains were used, perhaps the cradle motion could be un- 
coupled Srom the Lulu and the Alvin could dock with a virtually motionless 
cradle at 100' and then be peer to the surface rapidly and safely through 


the air-sea interface. 
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(ewe rete) Shailvsis of the Bulu-cradle interaction is the primary 
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To this time the ship cradle interaction motion has been studied using 
a simple model shown in Figure 2. While this model applies very well to a 
svstem with a single support it does not predict the motion of a four point 


suspenced cradle with ship pitch, roll, and heave. 
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Figure 2. 


Looking Yee Treure 3 it is apparent that the motions of a cradle be- 
neath a ship are extremely complex. The cradle can exhibit vertical motion 
(y), pitch (6), roll (4), fore and aft sway (a), transverse sway (8), twist- 
ing =bout the y axis, as well as horizontal motion in the x direction. In 
effect, all six degrees of freedom can be achieved with the cradle. For- 


tunately, to somewhat simplify the problem, fore and aft sway (a), 
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transverse sway (8), and twist were experimentally found to be insignifi- 
Sant . 

Consequently, the problem to ibe solved is a three degree of freedom 
one with roll (6), pitch (6), and heave (y), and where all three modes of 
motion can transfer energy back and forth thus contributing to the com- 
plexity of the problem. | 

Before beginning the modeling the initial conditions must be defined. 
To begin with, the cradle is assumed to be of length 1 = 17.5 feet, width 
w= 16.3 feet, with a weight of approximately 508 slugs. Since the 
reterence material has indicated that a possible mode of operation would 
be with a flat plate affixed to the bottom of the cradle, this is the mode 
chosen for our model. Consequently the eradie is assumed to be a rectan- 
gular flat plate with mass equivalent to that of the cradle alone. In the 
final analysis it will be shown that the added mass terms associated with 
the cradle so overpower the actual mass that the exact configuration of 
the cradle structure is inconsequential. aa 

It was first necessary to pick an adequate co-ordinate system which 
would accurately describe the system and make it possible to completely 
position the cradle knowing only three co-ordinates--heave (y), pitch (@), 
and roll ($). A care: system (Figure 4), fixed above the cradle was 
chosen with X positive forward, Y positive down and Z@ positive to port. 
Using the conventional right hand rule, positive 6 and positive ¢ are as 
shown in Figure 4. As is conventional in dynamic analysis the cradle is 
given an instantaneous dispilacenent positive y, positive 6, and positive 9. 

Since the equations of the system which are being dealt with are 


quite complex with many terms, it will be convenient to set up the 
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Cradle Equilibriun 


Figure 4, 
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resultant equation in standard state variable form consisting of six first 


order differential equations and an output expression. The most general 


form of such an equation is 


X,(t) ab Ajo eae Ayn X1(t) by by90--+-Din U,(t) 
X(t) iD... X9(t) : ; 

= , + : : 
X(t) Pens ee jee (58) Dayle ooo osm ORCS) 


oY 


X(t) = AX(t) + Bu(t) 


. 


where A is the N x N system matrix with B an N x r control matrix. 
strictiy speaking, this above representation applies only to a N- order 
linear equations, however, since there will be a multitude of terms, both 


linear and nonlinear, it will be a good bookkeeping aid to arrange terms 


as such. ., 


The basic equation of motion for a spring, mass, and damper system 


is used here. Development of the equations of motion is begun with the 


one 


consideration of this free body diagram. 
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where: x forces = ma WEIGHT Figure 5. 


or: 
my = - Fdrag - F spring + force (exciting) 
and is the case of three degrees of freedom 


MY = - Fanaa - fenning Ft feyciting LV 
ww —_ 


1)8 = - Parag 7 spring + Fexciting (85) 


Ino = - Farag - Fspring + Fexciting (9). 


Each term will be discussed independently. 


The Spring Force a 
The four springs (chains) at the corners of the cradle are the major 
factor in decoupling the surface motion from the cradle. The spring force 
is proportional to its stretch and this stretch may well be different for 
each of the four springs. This is one of the mechanisms by which energy 
is transferred from one mode, such as heave, to pitch or roll. Ideally, 


if all four springs were equal and the cradle were given a pure pitch, roll, 


or heave motion, then the other modes would not be excited and you would 
i 





=O 
have simply three uncoupled systems. However, considering the inexactness 
of the Peal World, the cross coupled damping, and the exterior forces, it 
would be unrealistic to perform such a simple analysis. In actuality the 
ship pitches, rolls, and heaves about an equilibrium position as does the 
cradle. The stretch of the four springs is defined as mre dtemerce ne 
position between the position of the ship and the position of the cradle 
relative to their respective equilibrium for each motion mode. 
The re erive Stretches of the four springs Ki> Kyo, Kg, and Ky are 


Ao) toOllows: 


6, = (y-Yosinw)t) + (L/28-L,/28 sinwt) + (W/2$-W /2%osinw st) 
84 = (y-Ypsinw,t) + (L/26-L/2@)sinw,t) - (W/2>-w/2e>sinwot) 
On = (y-Yp9Ssinw,t) ~ (,/29-L/28)sinw,t) = (W/29-W/22,SinwWot) 


6, = (y-Ypsina,t) - (L/28-L/26,sinwt) + (W/26-w/2%sinw,t) 


This is assuming that Yo,;6,, and Oy» are the ship's pitch, roll, and 


0? 
heave. 

Carrying ae a bit further, when the stretches are multiplied by 
their spring constants, we get the respective restoring spring forces and 
moments in general form. 


The equation for heave-- 


= i 6 


my = = Kp oe er 
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The equation for pitch-- 
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heave equation-- 


Plugging in: 


—— My eae W Nee gas) eas 

5h = KE Cyy - Yosinw,t) + & a 7 9)S8inw,t) + (> 5 Srey Sinw,t) J 
L : W Wile ew 

a Kol (yy - Ygsinw,t) + (3 Ys ~ , 898inw,t) - SiG - 7 %)Sinwt) J 
; L ii . W W , 

ays = Ygsinw,t) - (5 ¥3 Seas - (Ys - 5 PoSinwet) J 
e is Ties ees W ee a. 

oe Kyl(y, = Y)sinw,t) - oer - > SoSinw,t) + (> ys - 5 %Sinuwst) J 


pitch equation-- 
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roll equation-- 
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DBeavang Force 

The force which excites the cradle is transferred through the springs 
ferene cradle. In this analysis it is setasa that Lulu is heaving and 
pitching sinusoical with the same frequency w), and amplitudes Yo and 80> 
respectively. Additionally, it is assumed that Lulu is rolling sinusoidal- 


ly with amplitude 2 and frequency W». 


. Drag. Forces 
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terms to analyze since due to the 
low Reynold's numbers involved the drag forces are primarily caused by 
hydrodynamic form drag rather than viscous drag and are nonlinear. In ad- 
dition, for geometric shapes with hard edges the drag coefficients are 
generally independent of Reynold's numbers. 

The general equation describing drag forces on such a shape is given 
as Drag = eepcd V , and it 1s im this where we fimd@thesprobiem., If 
drag is to always oppose motion as it does here, the equation above should 
actually be written Drag = dep A.calv[V, which will preserve the correct 
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differential element of drag is defined in quadrant (1) of Figure 6 as) 
Bidras) = s0Cpdxdz (y+x8cosdcosé+z¢cosdcosé) |y+xdcosécosd+z$cos¢cos6 | . 
then the total drag on quadrant (1) would logically be 
% W/2 
D2AG = if f 4p Cpdxdz(y+xdcosdcos$+zdcos@cosd) |y+xbcosécos$+z$cos8cos9| . 


0 0 


“mis eprroach could be pursued to an exact solution if an identity for 





fe ba o 
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A+3+C] was obtainable. However, (A+B+C) | A+B+C| only equals 
(A+3+C) (A+B+C) when (A+B+C) is a positive number. The determination of 
whe exact sign of (A+B+C) or (y+x@cosOcos$+zdécosécos$) can only be obtsined 


fee & *S.Can Value of X, or Z is available. Consequently, if X and Z are 


+7. Tne integrand they cannot be assigned a specific value as they vary from 





ME et eee PeGEoCclively. Thus, as maihematically satisfying as @a 
Mace Cquction corm drag would be, it is alus, impossible to obtain. Due to 
this, a more straight forward yet more tedious finite element method was 
used to epproximate drag forces. 


The cradle was first divided up into four quadrants as shown below: 





Fbeniaeeen 7. 


Using the sign convention previously described, the velocities of any point 


in the four quadrants are as follows: 


v5 = (y + x 6 cos8cosd + zbcos@cosd) 


Vo = (y + x @ cos8cosd - zbcos8cos) 
V4 =(y - x 6 cos8cos¢ ~ zbcos§cosd) 
Va =.(y - x 9 cos@cos$¢ + zdcos8cosd). 


Each quadrant is then divided into sixteen areas and the drag is calculated 


for each one as shown for element 255 below. 
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Voo = (y+ = L6cos8cosét <= Wécos®coss) 





(Quadrant 1) 


aicume oy 


Wo. 


: coe cae : oa 
Drag(a5,) = 4 PCy (s a) (y+ jerecosecoso+ [E¢coshcosd) [yt 16 uecos8cos> 


ee 
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Hence, the velocity is assumed to be fixed over the entire element and thus, 
SO 1S the $18n OL tne arag. LaeaLlty, aS une tumvuer UL Ciemelius guES up, Luc 
possibility of the total velocity changing sign in an element goes down un- 
oe in the limit, you have numerically integrated the whole surface to an 
‘exact result. In this analysis 64 elements are taken each time step and 
added together to produce the nears drag. 

Pitch and roll drag are calculated in much the same way but with the 
added complication that ane drag element acts through a finite moment arm 
ee produce the drag in terms of foot lbs. One difference is the proper 
description of the velocity of the element which for our same element apo 
would be 


= : we) ° 3 th 
Vaoo = (ycos6cosd + oF: L 6+ ze WO). 


The computer subroutines for calculating the drag forces are given in 


the Appendix. 
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In cases where the relative velocity between a submerged body and its 
Ssurroundins water does not remain constant, additional forces come into play 
as a result of the acceleration of the mass of water entrained by the body. 
This added mass depends on the physical dimensions of the body as well as 
its mode of motion. The actual body mass plus the added mass is called the 
virtual, induced, or total mass. Similarly, there are added moments of 
inertia terms which accompany the body's angular acceleration. Conse- 
quently, for a long flat plate of dimensions such as 

|__| 

—OOO" eee 
the added mass of the entrained water for motion perpendicular to the plane 
memcmomcradie 1s definec as M, = oma*, or in our case a = 2 depth = L se 


Me = komW*L. The total mass is then actual mass plus added mass 


aE 
Mp 


The added roll moment of inertia due to the entrained water is defined as 


M(actual) + M (added) 


508 + yonw-L. 


or in our case Tap 1/89 1a" , 


eS 1/8o0(5) 41. 
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The total roll moment of inertia is then 
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Sammarly., for patch it is 
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These forces conzisting of buoyancy and weight are constant and thus 

@o not enter into the equation in any way. 
Overview 

The equations developed to this point are nonlinear due to the drag 
terms but fortunately can be solved by various numerical methods of inte- 
gration. The equations, arranged in state variable form are shown in 
Pigure 9. It is important to note that the drag terms shown are extremely 
complex and are functions of five of the six state variables. More specif- 
ically, the heave, pitch, and roll drag forces are functions of the follow- 
ing: y(2), vertical velocity of cradle; y(3), the cradle pitch angle; y(4), 


the cradle pitch angular velocity; y(5), the cradle roll angle; and y(6), 
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drags are different functions of these variables. 
The forms of these functions have been explained in previous para-~ 
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Figure 9, 
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DEVELOPEENT OF LIIEAR APPROXIMATIONS 


As discussed previously, the equations as developed are nonlinear. 
Mntortunately, ee femperesulscamin the necessity of a brute force ap-— 
proach to the solution. That is, the well developed test and analysis 
methods wnich apply to systems of linear equations ada be used to draw 
general conclusions about a nonlinear system. Consequently, linear ap- 
proximations for the drag terms would be helpful only if they were shown 
to give reasonable results when compared to the exact ioe” 

Ideally, the nonlinear drag terms, which now are functions of five 
of the six state variables, should be reduced to functions of only one 
velocity term. Actually, it would be reasonable to expect the majority 
BeDarLiCUler Urug veri LO ve caused vy Ule veauLi ly wean, Gl etashet 
With that term. That is, heave drag should mainly be caused by heave 
velocity, pitch drag by pitch angular veloc and similarly, roll drag 
by roll angular velocity. | 


Therefore, using this as an initial approach, the single variable 


heave drag equation might be: 


Drag (heave) = 5 pCpWLy? 
The equation for pitch drag would be somewhat more complicated and would 


be calculated as shown in Figure 10. 








Figure 10. 


a (drag) = 2x pCyWaX(X6) (XO) X 
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DRAG (pitch) = ae 


Similarly the equation for roll drag is 
PCpHWL*d? 
PeAG Gaedd)) = —— ae 
; 64 
As yet, even these equations are nonlinear due to their squared velocity terms. 
Several approximations to velocity squared damping have been used in the 
past. One of the commonly used electrical approximations consists of equat- 
ing the energy dissipated by a current squared resistor to the energy dissi- 
_ pated by a linear resistor when both are driven by a Sinusoidal current 
source. When this is done the I? term is in effect related to a constant 
multiplied — Single power of I. The constant which results is = tga 
where I, is the maximum current value and w, is the driving frequency. 
__Yherefore, if we consider the sea as a velocity squared damper, a rea- 
sonable approximation to heave drag would be 


8 2 
1, —— 





= OGe 
where W,¥o is the maximum vVel@eei cy OL y. i= Chis saeproach sis applied to 
@even and roll the following equations result: 


h  OOCpLVYqw 1y 
37 


Drag (heave) 


eCpLW8qw 6 


Drag (pitch) = on 

PCW" Loy ¥5¢ 
D aes = 
rag (roll) ae 


To check the validity of this linearization approach as well as the ac- 
curacy of the 64 element method, a comparison was made with values obtained 
from the exact drag eauations shown in Figure 10. Since these exact equa- 
tions are only valid for cases where single velocity terms are used, the com- 
parison was done for a case where heave drag is only a function of heave 
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drag is only a function of roll angular velocity. In other words, simple 
roll, pitch, and heave drag were tested independently. 

The results of the comparison shown in Figure 11 indicate clearly how 
the finite element method compares almost identically with the exact values 
@s computed with the equations of Figure 10. The one difference, that being 
that the linear approximation gave lower peak values as they appeared more 
Sinusoidal in shape was anticipated because we are equating energy dissipa- 
tion rather than peak values. 

Additionally, to assure that the one term approximation will give the 
desired drag force at the proper time, a comparison was done between actual 
and linear drag values which would be created during normal cradle oscilla- 


tions. Figures 12 through 14 compare drag values with Lulu exhibiting a 
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rte @eocdee ie25 ft. heave amplitude. The coxoeri- 


ee 
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meee second OO 
Seescno«~s that for this varticular driving frequency and amplitude the linear 
approximation is in phase with the actual values, yet tne amplitudes do not 
correspond. Another comparison at twice the frequency of the first confirmed 
that while a single variable drag term approximates the timing of the drag 
force well, it does not predict the amplitude. 

In an attempt to better estimate the drag forces several two and three 
term approximations were investigated but the results showed that, in gener- 
el, complete confusion reigned. Even worse, the drag forces were no longer 
in phase with tne actual values. 

The incorrect drag amplitudes result mainly from the terms in the coef- 
ficient of proportion of the linear drag term. In particular, for heave 


ty 


von, denends heavily on W.. the freauency of the cradle, 


Vee = = +7 cea em tant 


- 
WA od a) wre es ween See 


mae Ya_, the maximum displacement of the cradle. The cradle frequency art 
man De assumed to be the same as the driving frequency, but Yo, the maximum 
cradle amplitude, cannot be assumed to be ne same as the maximum driving 
Beeirtude. mm acttalaity, Yj or “bie cra@ie is aie maximum displacement of the 
Sredle as given by the equations of motion. However, that value of Yo9 de- 
pends on the chosen value for the @racg terms whieh in turn is@a funct#on of 
Yo: Pesiically wewNaverdn 1teracave problem where values for the drag terms 
must be initially assumed and comparisons must be made with the actual non- 
linear responses. 

To provide a baseline for comparison a brute force analysis was con- 


ducted on the nonlinear eguations to determine the actual response or gain 


of the system at various input frecuencies. Ten separate frequencies were 


analyzed and the steady state gains were plotted versus frequency. The same 





es 
-30- 
Mumeols was re.eaved but this time using the single term linear drag an- 


proxinations in place of the accepted 64 element drag values. This analysis 


wes t1so ccnducted at the ten different driving frequencies with three dif- 


% 


ferent driving amplitudes. 
Figures 15 and 16 show clearly the differences between the linear and 

nonlinear responses at various driving frequencies. The roll response curve 

me Ommitted since it falis just Serge O11 beins On UCpmem Che plticn curves. 

These figures show that the linear approximations peak at the same frequen- 

cies as do the actual response curves but the gains vary widely depending 

on the input amplitude. In general, the linearized response curves show a 

condition of lesser damping then is actually the case. However, as the iis 

‘ing amplitudes go up, the apparent damping increases due to the presence of 
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SEALS ws reveyyv 


He 
© 


the larger value of driving amplitude (Yj; in the coeff 
-Since the main value of such an analysis is the prediction of the 
resonance peaks, and since the linear and nonlinear responses are almost 
identical except in the eA eeOus . Lee 250: — alae is valuable to correct the 
n 

linear response curves by increasing the apparent damping until the linear 
responses conform to the nonlinear responses. This can be done by multiply- 
ing the previously discussed linearized drag terms by a constant depending 
on the magnitude of the driving amplitude . 

Comparing the linear response curves to the simple single degree of 
freedom system shown in Figure 17, we rue that the gain of such a system 
at resonance is given by the equation E sinere - is the damping factor. 
Using this equation we find that the nonlinear heave resonance peak compares 


well to a linear system with a damping factor of 0.193. Similarly the pitch 


response compares well with a damping factor of 0.167. After finding the 
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Figure 15. 
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Sinele degree of freedom forced vibration response, 


Figure 17. 
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indicated linear values of demping factors for the linear response curves 
With various input amplitudes, a plot of drag multipliers versus amplitude 
input was made. This plot, shown 7 Figure 18, eee the factor by 
which the dreg terms must be multiplied in order to make the linear response 
correspond closely to the nonlinear ene for each input amplitude. For 
example, for our normal 4.25 ft. heave and 10° roll and pitch amplitude, the 
heave drag term would need to be multiplied by 2.7, the pitch drag term by 
yeh, and the roll drag term by 4.0. 

Summing up, linearization of the nonlinear drag terms is a complex mat- 
ter. As in all problems in the physical sciences it is never feasible to 


express such linearizations with absolute certainty and with the inclusion 


~ 
_ 


of all possible variables. In fact, to treat any physical problem in a 
Su@reteioat Way io iS Mecessar, TO first make acsumptions about the system's 
behavior that may not be entirely correct. In this case drag is assumed to 
be linearly proportional to velocity, but, as is later discovered, it is also 
related to the driving amplitude and therefore appropriate correction factors 
must be applied. In spite of the difficulty involved, the value of this 
linearization is to provide the coefficients for the linear A matrix which 


can then be operated on to find the general trends of the system. 
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Figure 18. 
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Cineiik IV. 


SO morrOl@Or EQUATIONS 


The complete solution of differential equations may be split into two 
parts. a steady-state solution and a transient solution which dies out at 
a rate depending on the degree of damping, mass, and time. Our main interest 


here is with the steady-state solution. 


Time Domain Analvsis 

The time domain solution of the equations of motion was found by the 
use of numerical integration. The integration was done using the newly in- 
Stalled Interdata computer facility at M.I.T.'s Me chanical Engineering De- 


~ 


partment. The same department's program DYSYS, (DYnamic SYstem Simulation), 


, - ~ oi 8 “ ° oe ane Ve Se I ee oS 3 
moomoo d L£Or te Cille Ullal aiarvsis oi Ties FOR El Tete Ee: SHS bycvs 


Mives a set of first order linear or nonlinear differential equations and 
can be particularly useful for obtaining transient as well as steady-state 
solutions to nonlinear time varying problems. DYSYS uses a fourth order 
Runge-Kutta integration technique to obtain an approximate solution to this 
equation-- 


Given Yo> the value of the Y vector at some value of t, say ae then 


the approximate value of the y vector at ty + At is 


At : 
Bet atk * 2K t 2m Ka), 


the K's being equal to 


K = F(to + SS, ¥ + +3 K 


1 0 Qo) 





ee 
ae At ae 
en a 5 Yg + “K)) 


K. = F(to tees foe t KS) 


0 


The Runge-Kutta !‘ethod is derived from a Taylor Series Expansion about the 
point ty and Yo: Theoretically, thewmethod will converge to an exact solu- 
tion as t approaches zero.. The methcd has been shown to be stable, self- 
starting end to have a built-in error ee 

Since the program is set up for simultareous first order differential 
equations, the three second-order equations must be expressed as six first- 
order equations. This can be done by letting 

y =ey(i) == 


youl = y 


y(2) = y 
Similarly, 

OL) Sy 

aay, 

y(3) = @ 

y(4) = Q 

y(S) = 9 
-y(6) = 6, 


and hence y(1) through y(6) are the six state variables which completely 
define the system. 

The computer program which was written to effect the solution is given 
in the appendix. The program has as subroutines EQSYM, which indicates 


"Equation Simulation” as well as three separate drag subroutines which must 





=) 
be called each time step and calculates the drag at each point in the 
program. A time interval of 0.2 seconds was selected as a time step and 


Bave reasonable results. 


Frequency Domain Analysis 

Since not only specific solutions for given parameters were desired, 
a frequency domain analysis resulting in Bode plots was undertaken. Bode 
plots clearly illustrate the relative Folic of a system, with gain and 
phase margins often being defined in terms of Bode plots. 

This analysis was conducted using the M.I.T. Mechanical Engineering 
Department computer program MBODE. This program calculates frequency do- 
main gain anc phase plots for linear dynamic systems (with constant coef- 


ficients) expressed in matrix form such as 


- x = Ax + Bu 


y= 62 pu 
where 
| A = the system A matrix 
B = the input matrix 
O = oe output matrix (usually diagonal ones or 


appropriate weighting factors) ‘ 


D = the direct coupled output matrix (usually zero). 


The fact that MBODE only operates on linear constant coefficient equa- 
tions is the reason we were forced to develop linear approximations for the 


nonlinear drag terms discusSed.in section III. The Bode plots which 
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Smited from this method of analysis show clearly the properties of the 


ystem and thes are presented in section V. 
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CHAPTER V. 


DISCUSSION OF RESULTS 


Time Domain Response 

For this analysis numerous trials were run at various damping levels, 
input frequencies, input amplitudes, and spring constants. The majority 
or the trials were run at the baseline conditions of 4.25 ft. ship heave, 
10° moll and pitch, and 20 “o> spring constant with a 10% deviation among 
the four springs. The 10% deviation-is the method by which the nonideal 
Mature of the springs was modeled and is the vehicle for energy transfer 
between the different motion modes. 

In a discussion with J. Mavor of W.H.O.I. it was learned that the pre- 
viously mentioned amplitudes with a 5-6 second period could be considered 
Sueece) shipimetienc in e stete fowr seer Figure 19 shows this realistic 
Situation of heave, pitch and roll and indicates that the transmitted ship 


heave is reduced from 4.25 ft. to about 5 inches while the roll and pitch 


are reduced from 10° to 2°. A total effective spring constant of 80-100 
Rpt 


rt is required to accomplish this decoupling; however, the analysis il- 
lustrates the fact that ship-cradle motion can be uncoupled if soft springs 
are used. | 

In Figure 20 the pitch and roll amplitude have been doubled and the 
ship has a slow 10-12 second period to give a "worst case" example of pre- 
dicted cradle response. Although this decoupling is not as good as with 
the shorter period, the heave has still been reduced to less than one foot 
while the roll and pitch have been reduced from 20° to 10° and ‘im respec- 


tively. 
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Mrequency Domain Resnonse 

Since specific response curves do not reveal much about the general 
nature of the system, a frequency domain analysis was conducted in order to 
get more descriptive information on how the system will react to various 
input frequencies. Figure 21 is a plot of the system's gain response at 
baseline conditions of amplitude and spring constant but at ten different 
Pnput ences, It is, in effect, a Bode plot of the nonlinear case. 
The plot shows clearly the three reasonances of the system, the leftmost 
peak being the heave resonance and the rightmost peaks being pitch ad neon ell 
resonances, the roll resonance being the smaller of the two. 

mais figure required ten separate solutions of the nonlinear ‘equa- 
tions of motion involving great expenditures of Sour esas. Although 
the piot shows the response of the system under fixed conditions of spring 
estan, mass, and driving amplitude, it is clearly not a desirable way in 
which to find how the system responds when differing sea conditions exist 
or when the system parameters are modified. It is here that the value of 
the Mieevization process comes into full view. 

In Chapter III it was noted that the nonlinear see terms could be 
linearized once the driving amplitude was known. Again using the 4.25 ft. 
heave with 10° roll and pitch as the hegeiltsas we can find the linear drag 
approximation by referring to Figure 18, Chapter III, and Shuster the 
proper multiplier. Once these drag coefficients are set, Bode plots may 
be initiated which show the response of all outputs versus all inputs over 
@ wide frequency range. Additionally, the effects of variations in’ spre 
constant or mass may be sinshed With only minor enter ions to the now 


linearized A matrix. 
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Figures 22 through 25 show the system's response across a range of 
mae to 1.0 Hertz or a roll pitch and heave period of 1.0 to 100.0 seconds. 
The outcuts for heave versus ship pitch and roll show that the ship's 
mecching and rolling couples energy into heave motion which then peaks at 
heave netural frequency as well as the pitch and roll natural frequency. 
This fact indicates the dangers of designing a decoupling system to de- 
couple only heave. It is clear from the plot that a decoupling system de- 
Signed only for heave may well excite the pitch and roll modes at their 
natural frequencies. This fact was overlooked by single degree of freedom 
analysis. | | 

Several other conclusions may be drawn from Figures 22 through 25. 
While a 20 1bf/ft was chosen as a baseline for the spring constants and 


” e bn | 


igures ce and 23, it would be desirahie to predict eusily 


WaS Used ror 
what would happen with higher or lower spring constants. Figures 24 and 
e) show the response of the system with a spring constant of four -times 
that of Figures 22 and 23. It is not surprising to see that the resonant 
BP ixencies are shifted upward by a factor of approximately two. This in- 
dicates that in spite of the nonlinearity of the parent system, the system 
obeys the basic laws for calculating natural frequencies of translational 
and rotational systems. Being more specific, heave natural frequency is 


proportional to ——— , while pitch and roll natural frequencies are 


proportional to vKL and vKL respectively. 
ey: Ie 


As a result of the above, we have information on the system's behavior 
over a wide range of frequencies without the need of calculating a Bode plot 


for each value of K. In effect, using K = 20 lbf/ft as the starting point, 





=S0— 


we need only shift the response peaks up or down the frequency scale by 


a Wriew 
a factor of the square root of Kew! Koo» or c= | Wola AS 
K20 


This same method of peak shifting can be applied to variations of mass or 


moment of inertia in a like manner, but in this case the shift will be 


YMt(base) 
proportional to a or pauicoe!) 


summing up, this analysis has shown that motions of a submerged cradle 
can, in fact, be nearly decoupled from surface ship motions if a soft ef- 
fective spring constant is provided at the four cradle corners. In addi- 
tion, the analysis has shown that the critical design consideration is the 
pitch and roll natural frequency rather than the taal frequency of 
heeve. Ac = resguit. careful thought nee pe given to the sélection ¢c 
Spring constants since the equations of motion are coupled and are capable 
of transferring energy from one mode to another and hence exciting other 


modes at their natural frequencies. 
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CHAPTER VI. 
Pee CATIONS POSSIBILITIES 


fms NOt the Major Objective of this thesis to offer a detailed de- 
sign of a decoupling system. It is enlightening, however, to discuss 
several approaches that could be applied in the design of such a system. 
The primary consideration vere is to obtain a soft spring constant and 
Still support a large static load. 

From the previous analysis it is apparent that the soft spring con- 
stant must be at the four cradle corners. Further review indicates a 


spring constant between 20 and 50 1b,/ft would seem to be a reasonable 


~ 
- 
ee, 


selection for these cormer springs. 

ft thic voint there are two tvves of approaches to conmsider. The 
first and simplest approach is the passive system where we try to create 
a spring with the proper softness so that without any outside control the 
Springs will decouple the ship-cradle motions. This approach can be 
likened to vibration isolation of a motor or a vibrating mass. It is sim- 
ple and reasonably foolproof. 

The alternative method is of course a more complex active system 
where some sort of control is used to actively uncouple the two motions. 
Mis control can be in tie form of position, velocity, acceleration, inte- 
gral or differential feedback, horizontal sensing transducers or mechani- 
cal linkages. Consequently, unlike passive systems, these types of sys- 


tems can be tuned to changing sea conditions or conditions of loading, 


thereby giving them a possible advantage. 
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passive Sys ems 

To create a passive system for the cradle would require consideration 
of several methods of achieving the desired soft spring constant. A first 
choice for component selection would probably be simple metallic springs 
which deform in accordance with Hooke's Law. Unfortunately, the static 
cradle weight of about 16,000 pounds would stretch the four 50 ie ot 
springs to about 80 feet before reaching equilibrium. Clearly, this would 
be 1 SERCO for if coil springs were used they would necessarily re- 
semble long "slinkys’' and would prove completely unmanageable. Leaf springs 
would similarly be of little use as they too would deflect excessively in 
the static mode. 

Torsional springs would possibly work if they were rigged through a 
complex gearing system which wound up as maa was abpiied; howéver, they 
too would require a large static deflection under the simple cradle load 
and would involve complex machinery and cable drums. Such a mechanical sys- 
tem would probably be more of a liability than an asset in an ocean environ- 
ment. 

Another possible alternative might be to take advantage of the non- 
linear stress strain properties of elastomers. Vulcanized rubber will 
stretch to nearly ten times its own length and support up to ten tons per 
square inch. Unfortunately, the resistance to deformation usually increases 
in greater oroportion than does the extension. If, by special treatment, an 
elastomer with the opposite properties were developed, it might be a logical 
choice to use as a spring. However, a short literature search turned up no 


Suggestion that such a material exists. 
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A more reasonable method of producing such an effective spring constant 
feuld be with a hydraulic or pneumatic device. A device such as a hydraulic 
or pneumatic cylinder has the advantage that it can be biased to an initial 
static load and then oscillate about that equilibrium position with an ef- 
fective spring constant which depends on pressure, temperature, properties 
ee the fluid, and dimensions of the system. 


To illustrate such a System consider Figure 26. 


| FORCE 
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Figure 26. 


This system would have a spring constant which would depend heavily on kK, 
the compressibility of the fluid. It is evident that a fluid with a high 
compressibility, such as air, would result in a softer overall spring con- 
Stant which is desirable in this case. Additionally, a decoupling system 
with alr springs would be capable of being set to an initial static force 
equal to that of the cradle weight. This might be accomplished with a sim- 
ple pressure regulating device. Following initial set, the cradle could 
‘then oscillate about equilibrium with a soft effective spring constant de= 


pending on the previously mentioned parameters. 
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Before considering this simple approach there are two important con- 
Sideraticns to be aware of. First, the oscillating piston must look exact- 
imeerike a soft spring to the disturbing ship motion. That is, the differ- 
ence in piston force from equilibrium to the extremes of piston travel must 
compare well with what would be seen if a real metallic spring were being 
msed. in other words, for a ten foot pneumatic cylinder with a normally 
centered piston (capable of + 5 foot excursions), the fully extended piston 
force nee be equal to static cradle weight plus K (effective) ° 5 feet. 
For example, if the desired spring constant were 50 Mor) Ets then the piston 
could only exert 250 lb, more than the cradle weight at its maximum stroke 
of five feet. Quickly calculating the change in pressure for full compres- 
Sion using P1Vj = PoV5, it is obvious that if the pneumatic system were 
gead endea, the reduction in vuoiume to 1/5 of itS preVious Valué would re- 
sult in a fivefold increase in pressure. 

Consequently, the first consideration in such a design is the fact that 
@ plenum chamber would have to be used in conjunction with such a system in 
order to nae the desired results. To illustrate the size of the required 
plenum chamber, we can again use PiVy = PoVo (if se To). Pecunia a base- 
line piston area of 1 ft.2, a static cradle load of 000 1d > per spring 
would require a piston pressure difference of 278 psi in order to maintain 
the cradle at equilibrium. Calculating the plenum size required with the 


following conditions 


AVeat fuliostroke = oert.- 
AF at full stroke = 250 1b, 
P1V) = PoVn or 28 Va emiOn te tony ~ Seca) 


3 
V, = 815 ft 





=b1= 
me uniorturately find it is a chamber far exceeding the capacity of the 
Bracie to easily carry. 
However, further calculations show that the size of the required plenum 


Will go down as working pressure goes up. For example, 


assume PL = 2000 psi 
Piston area = 2 in.? 
or POCCm ie 2000 +125). (V. =. .069) 


a 
= 5 
Vy act 9.4 feta 


a 


and is a volume which is well within the carrying capacity of the cradle. 
The second major consideration in passive system design is the problem 
of centering the system. With a pneumatic spring such as is sketched in 


Figure 27, the piston and plenum might be pressurized on the surface until 


the piston was positioned at midstroke and supporting the cradie's weigh 


PISTON 





va 





sae OD, 


Preure 27. 





aho> 
Mec einivial Setving, @lthough correct on the surface, would be incorrect 
after the effects of hydrostatic pressure, temperature, and leakage had taken 
place. This fact would necessitate the constant bleeding in of make-up gas. 
In effect, some active control would be necessary to make this passive design 
workable. 

In addition to the aforementioned major design considerations, when us-— 
ing a gas, with its higher compressibility, there are serious technical prob- 
lems to be overcome. These problems involve containment of a gas under high 
pressure, sealing, and danger of explosion, all of which may lead to bulky 
and expensive design. 

Throughout this entire first section only passive decoupling systems are 
discussed. These systems all have the same characteristics in that they os- 


bo] 


a i ¢ eect) ee sea os maaan ts Ce ad = eat = . Vee — ~~ oe om =} 
ciliate ahnout an ediiilibrium position not determined actively by sensors and 


q 


feedback, but by cradle weight, system dimensions, and initial conditions. 


Active Systems 


oystems which employ active feedback for control are called active sys- 
tems. Typical examples of these systems include automobile power steering, 
hydraulic positioning devices, and automated milling machine tables. These 
meiive systems are all characterized by a mechanical, electrical, or hy- 
draulie system by which information on the output is fed back to control the 
actions of the input. This feedback control can be designed to provide fast 
or slow response, maximum or minimum overshoot or undershoot, and nearly any 
natural frequency within the constraints imposed by the parameters of the 
system. 


The simplest form of an electrohydraulic position control device is 





Soo. 


Bhown below in Figure 28. 
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Figure 28. 


The feedback Loop cons1Sts OI sensors wnicn leeu vacn pusiviuu, vesuc- 


ity, and acceleration information which improve the system's response to 


disturbances at the input. Position feedback often is sufficient for system 


control but does not provide the tight feedback loop that results when 


velocity and acceleration are also used. 


A typical response curve for such a system is shown in Figure 29. 
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Figure 29. 
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mhaet if, the system would tend to rise to the desired level, overshoot it, 
and then oscillate about that level with a decaying oscillation edie: 
on the parameters of the system. It is this oscillation which might well be 


used to create a virtual soft spring constant at the constant load. 


A Proportional Control Approach 


As previously stated, feeding back all three variables would tend to 
make the system have less overshoot and en ce Since we are attempt- 
ing to utilize the system's natural oscillations, we would be more success~ 
ful using only a position feedback system which, while not enhancing oscil- 
lations, does not discourage them strongly. 

To determine the feasibility of such a system and to get some design 
Meencters for guidance, the system equations are determined below. 


First, looking at the spool valve and assuming the hydraulic fluid is 


Water, 
Valve flow-- 
o> ‘Coe X - KopF'm 
were QF = Tlow to load 
Kox = flow sensitivity or valve gain 
Kop = negative slope of valve characteristic curve at 


constant value of x. See Figure 30. 


Figure 30. 
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Fer an ideal valve the slove, Kop? 
puch @ valve is called a flow control valve. 
Valve forces-- 
—-Kpx + I 0) 
where Ep = valve control spring 
PF; = valve force 


Next, considering the piston, 


Piston flow-- 
ee aay 
Qn = AY + Py 4s (Ke + gi) + CoP. 
where: 
Co = piston leakage 
ie cemplilance Of piston 
V = volume of piston chamber 
om = bulk modulus of fluid 


Piston forces-- 


ae & 
MY PA Wy 
A = area of piston 
M = mass of load 
W. = weight of static load 


Combining all four equations, 


ee V 
= 1, a 
oe ENE AON or) ail O25 





m 
and: My = Pat fas 
mimally , 7 
he nies C M.: 
_ V aun 
A Cre )) At Ay = ROX 


iS very small and can be assumed as zero.. 
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Mec using Larlace transforms, 


oe) 
y(s) = Cai 
» MI Ja V o 
Ss ( S ft Ne 2 z,” + S Ae te A) 


such a er can be described in terms of poles and zeros and valuable 
information can be obtained with an S plane plot of the position of these 
poles and zeros. Our case results in three poles, one ne the origin and a 
pair of complex conjugate aren The conjugate pole pair indicates the 
system is oscillatory while the pole at the origin indicates that the system 


may not reach the desired new setpoint or that it will have a steady state 


Deere seowewer, this 1S Of dittle concern to us here. Plotting the zeros 


im Figure 31, YZ \ 
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Figure 31. 


The natural frequency of interest is defined by the length of the 
phazor connecting the origin to a complex pole. Basically, if we put the 
transfer function in the form 


T = const 
m7 s(s + a, + J8,)(s + a - §J8,) 


3 


we can adjust the natural frequency of the system by adjusting the coef- 


ficients which control the position of the complex conjugate poles. 
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In our case, 


M V 


calling 5 Ke + a Ao 
M _ 
8 
A = Ay 
and rewriting 
m const 
y(s) 


s(A,s° + Ays - Ag) 
The natural frequency of the system is defined by the roots of the 


equation Ay Ao 


or the megnitude of 


or Al 2 A? 
Py + GY - wn 
h 


= 0 and all terms with A 


If we assume leakage is equal to zero then Cy 1 


go to zero, and the expression reduces to 


—-4uA 
a 
n m V 
boa Ke + —)) 
D 
. OA 


n MCK, + a 


To get a natural frequency of 0.03HZ or about 0.2 the argument of the 
square root would have to be about 0.50. 
Assuming a baseline piston size of 1 ft? and V = 10 ft? using a mass of 


2000 slugs and the bulk modulus of water B= 320,000 psi, the equation gives 
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1 2 ee 
; Ze sone felt LO 
0.90 Secé 2000 reise (kK, + 
52,000 
er. With vroper units 
— (144)(144)(2) int 
ES sect USF) = See site 17280 in? 
DOO) a ee ee 
Tt 12in 320 ,000psi 
in? 
Ke =. 00 —— 
n° 
Hence a compliance of about 500 psi must be built into the system. Assuming 


an air accumulator with a working pressure of 278 psi and change of one psi 
ey ee Iti 
Pie = 20 (V3 - 200) 
V, = S0uft> i 
whieh acain ig verv jiarge and noikv. 
Fortunately, as was shown in the passive systems, this large accumulator 
may be made smaller by use of a higher working pressure. For example, 
assume Py = 2000 psi 
Area = 2 in? 
a) in? 
0.50 = 


2000 slugs lft 
22 in 


20 
(Ke + 359-000 ) 
» 3 

fh 10,006 = 
& Dol 


Or again, using an air accumulator at 2000 psi, one psi should change the 


volume 3 
in 
psi 


0.096 


PV, = PoVo 





e510 


2000 V 


1 eOoiey 


‘a 0.096) 


192 in? 


tH 
<4 


which is of reasonable size and is easy to build. 

Actually, a working pressure somewhere between the limits discussed 
here would be most appropriate in order to reduce piston size and weight 
as well as reduce the plenum size. In short, there is probably an optimum 
selection of working pressure, piston area, and hydraulic fluid. A no- 


tional system employing this approach is shown below in Figure 32. 
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Pyetmcreasimemr@eeutevye)] Of comerol is include integral as well as 
proportional control we may be able to improve the system and possibly re- 
fimiee the Size of the plenum chamber. 

Modeling the system as shown below the following equations result, 


Figure 33. 
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INPUT COD 


Pp. = Load pressure 


9 | Kee — loval tiurd pilus 
' cylinder compliance 
an WLLL: 

| ee Q, = Flow to load 


Fosinuwt = Disturbance 
2 


— DRrannvtstanxnal ia 
R = fe Owos AS ade a con=- 


constant 


S = Integral control 
constant 


D = Damping constant 


ae Ket m a 


(2) My = -P,A-Fosinut-Dy 





(3) Q, = Ry+Sfydt-c(t). 
Combining equations (1), (2), and (3) 


My = - P - Fowsinwt - Dy 


but . ’ y Qn 


=o +} 


-Plugging in 


ose y Qn s 
i Ae as Met - Fowcoswt - Dy 











=e 


BESO, 
Q, = Ry + Sy ~ C(t) 
or finally, a ” ee. oe Sy : ~~ 
Peete yet > y+ = C(t)+Fow‘sinut. 
Using LaPlace transforms, 
<1 5) C(t) + Fow*sinut 
ee ae RS 
(Ms' + ps3 + 524 + _S 
Ney Key Ket 


Where C(t) is the input setpoint (time rate of change) and Fow*sinwt is the 


disturbance, the characteristic equation is then 


A S 
Hose soe ee = 0 
or ; nag A , me ar a. | 


uw) ~ 








St + — + oe oo + 
M Ke+M KetM KetM 


This equation can then be analyzed with standard procedures and proper co~ 
eiticrents for X, 5, end K.; can be found. 

some general comments can be made about proportional-integral versus 
simple proportional control. Addition of integral control allows the sys- 
tem to respond to 2a change in set point with zero steady state error while 
it does not change the basic natural frequency of the ere Moreover, 
this control element provides an additional independent constant which can 
be used to more accurately tune the system for maximum decoupling effect. 
The availability of this additional constant may well offset the complexity 
added with an integral feedback unit. 

One fact to be aware of with the P.I. control discussed above is that 


the cradle has been shown to exhibit less than critical damping. Hence, 





Lee 
Pale wcharacteristic equation would probably become unstable upon the ap- 
plication of a reasonable gain and additional damping may be needed. This 
damping could be provided either by mechanical means or by additional feed- 
back such as velocity and acceleration. 
ihe Geverminagzvion of the optimal design for such a feedback control 


System is left for future work. 


Other Systems 


Basically, all other systems are combinations of the passive and active 
Systems. Some system modifications which may be desirable include: 

(1) Inclusion of horizontal sensing valves which tend to channel fluid 
to the low sides of the cradle thus pire norizontal feedback to an active 
Ee eein. 

We) scuusuloution Of rotary air or fluid motors vith gear racks tO 
Beevide the vertical position control. 

(3) The utilization of an active-passive system which sesteirs a posi- 


tion control to center a passive piston or rack. 
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CHAPTER ViIT. 


CONCLUSIC?:S AND RECOLAIPNDATIONUS 


As a result of the investigations carried out in this thesis, several 
conclusions have ont arrived at. 

1. <A four point passive decoupling scheme can substantially decrease 
ship heave, pitch, and roll transmission to a submerged cradle. 

2. The effective spring constants required for decoupling fall in the 
range of 20-50 lb,/ft. 

3. Each mode of oscillation will couple energy into the other two 
modes, hence the design of a decoupling system must take each natural fre- 
quency into consideration. - 

4h. Cradle hydrodynamic drag is a complex function of frequency, ampli- 
fude, and velocity and must. be cnosen to suit Che inuiviud@k oL VW UG wid 
is to be modeled. 

5. Conventional springs such as coil or leaf springs which will have a 
soft spring constant and also a large initial static load do not exist. 

6. Passive pneumatic springs, while meeting desired requirements, would 
require large plenum chambers in order to operate at safe and reasonable 
pressure levels. Additionally, these passive devices would be dif PareulGeaeo 
center in the cylinder. 

ion nh logteciwapproach for creation of Soft springs witnweinitial static 
eee would be an active feedback control system. A system, such as one 
utilizing proportional and integral feedback, could substantially decouple 
ship motions from the cradle, yet maintain the resultant piston oscillations 


geeub the center of the cylinder or set point. 
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In summary, both the theoretical and practical aspects 
of a cradle mounted decoupling system have been investigated. 
The value of this analysis is that it brings us to the point 
where we must decide whether the possible improvements in op- 
erating conditions warrant the considerable time and expense 
involved in making such a system come true. If the analysis 
had shown that a simple passive system would work and was 
readily available, further tests and construction would be 
easy to justify. On the other hand, since an active system 
seems to be necessary, extreme caution should be exercised 
prior to construction to assure that; first, proper feedback 
coefficients are found; and second, the necessary control com- 
Menents are available and are readily adaptable to the ocean 
environment. In addition to the above, care must be taken to 
insure that any resultant gat aimes design satisfys all ee 
ty requirements which may be imposed by WHOI, the Navy, or 
any other regulatory bodys. 

Other means of decoupling ship-cradle motion have been, 
and are being evaluated. The author hopes that this work will 
enable WHOI to more clearly evaluate the cradle aeaneee system 


and how it compares to other methods being considered. 
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